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Abstract. We review the formation scenarios for binary black holes, and show that their co- 
alescence rate depends very strongly on the outcome of the second mass transfer. However, 
the observations of IC10 X-l, an binary with a massive black hole accreting from a Wolf- 
Rayet star proves that this mass transfer can be stable. We analyze the future evolution of 
IC10 X-l and show that it is very likely to form a binary black hole system merging in a 
few Gyrs. We estimate the coalescence rate density of such systems to be 0.06 Mpc~ 3 Myr~' , 
and the detection rate for the current LIGO/VIRGO of 0.69 yr" 1 , a much higher value than 
the expected double neutron star rate. Thus the first detection of a gravitational wave source 
is likely to be a coalescence of a binary black hole. 
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1. Introduction 

Detecting binary black holes (BBH) is ex- 
tremely difficult. They do not emit significant 
electromagnetic radiation, like the black holes 
(BH) in accreting systems. They can only be 
discovered via their gravitational influence: ei- 
ther through microlensing or through obser- 
vation of burst of gravitational waves when 
they coalesce. Microlensing caused by a BBH 
is an extremely rare event, we know only a 
few cand idates for microlensing through sin- 
gle BHs dBennett et alj|2002t iMao et alj|2002t 
iPoindexter et al.ll2.005b . and the the microlens- 
ing by BBH must be much less frequent. The 



development of gravitational wave observato- 
ries on the other hand is very promising and 
may lead to detection of a BBH coalescence. 
The current sensitivity of LIGO and VIRGO 
interferometers permits to detect a coalescence 
of a system of tw o 10 M H BHs at a d istance ex- 
ceeding lOOMpc (lAbbott et al.l2009l) . The sen- 
sitivity of the instruments is being improved 
which brings us closer to a possibility of de- 
tection. 

In this paper we investigate discuss the pos- 
sibility of existence of BBHs. In section 2 we 
outline the challenge of formation of a BBH 
from the point of view of binary evolution. In 
section 3 we summarize the observations of 
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IC10 X-l, and discuss its future evolution. We 
estimate the merger and detection rate of BBHs 
based on IC 10 X-l in section 4, and in section 
5 we summarize and discuss these results. 



2. Challenges in formation of binary 
black holes 

The formation scenario of BBHs has been 
discussed in the literature by many au- 
thors, see e.g. Lipunov et al .1 (ll997alTbh; 
Belc zvnski et all (120001); 



Grishchuk et al.l 



d2001l) : lBelczvnski et alj d2007D . Here we pro- 
vide an outline based on the StarTrack binary 
popul ation synthesis code (iBelczvnski et al.l 
120021 120081) . 

In order to form a tight BBH with a merger 
time shorter than the Hubble time we need 
to start with binary consisting of two massive 
stars, with the mass above 20 M . They will 
evolve very quickly and soon one of them, 
the initially more massive will enter the giant 
phase. This will lead to a mass transfer that will 
initially be meta stable but once the mass ration 
becomes close to unity it will stabilize. After 
this initial mass transfer the system will con- 
sist of a He star - the core of the initially more 
massive one, and the rejuvenated companion, 
that has gained a significant amount of mass 
from the companion. The Helium star evolves 
further and soon a first BH is formed. Now 
the system consist of a BH and the companion 
that will soon become a giant an start trans- 
ferring mass onto the BH. This ensuing mass 
transfer is the bottleneck of the BBH forma- 
tion. The BH has a mass that is typically much 
smaller than the rejuvenated donor. Therefore 
the mass transfer is not stable and the system 
enters a common envelope (CE) phase. The 
system may survive it yet if the donor has de- 
veloped a well defined structure of a core and 
envelope dTaam & Sandquisll20 00). The simu- 
lations, however, show that the donors do not 
have that structure and the CE phase leads to 
a merger rather than formation of a tight bi- 
nary. However, if a binary manages to survive 
this mass transfer episode than its further evo- 
lution is simple. The donor is stripped of its 
envelope and it becomes a He star. This star 
evolves also quickly and explodes forming a 



second BH. Thus if a system avoids a merger 
during the second mass transfer than a merging 
BBH can be formed. 

These considerat i ons h ave been presented 
by IBelczvnski et al.1 (l2007h . That paper con- 
cludes that although BBH are very bright and 
detectable by gravitational wave interferome- 
ters up to very large distances the effective 
merger rate will be quite low since the forma- 
tion rate of BBH is extremely small if the sys- 
tem do not survive the CE episode. However, 
our knowledge of the physics involved in CE 
evolution is poor and this leads to huge un- 
certainty in the calculations of formation and 
merger rates od BBHs. This conclusion was 
based on the analysis of binary evolution of 
stars with the solar metallicity. The evolu- 
tion of metal free star h as been considered by 
IBelczvnski et al.ld2004 . 

3. IC10X-1 

IC10 has b een discovered 120 years ago 
dSwifll Il889f) . It is an irregular d warf galaxy 
at distance approximately 660kpc dSaha et all 
119961: ISakai et al.ll 19991: feorissova et al.ll2000l) 
It has a rather low metallicity of 0.1 5Z 
(|Lequeux et al.llT9~79T) . and it is the only star- 
burst galaxy in the Local Group, with the 
star formation rate of .04 - 0.08 M G yr -1 
dHunter & G allagherll 19861) IC10 hosts a very 
large number of WR star s dRicher et al.ll2Q0lb 
iMassev & Holmes! |2002|) . with the surface 
density of WR star being four times larger 
than anywhere else in the Local Group even 
for the regions of highest star formation 
dMassev & Johnson! 1 19981: iRoyer et all 1200 it 
ICrowfheret alJl2003h. The R OSAT X-ray ob- 
servations (Brandt et aill997t) revealed several 
sources, with IC10 X-l being the brightest. 

3.1. Observations and the present state 

IC10 X-l has already been reported to show 
some X-ray variability in the discovery pa- 
per dBrandt et al.l Il997l) . The optical com- 
panion h as been identified as a WR star 
MAC92 (IClark & CrowtheJ [2004). with the 
use of WR star catalogue s obtained pre- 
viously (Masse yet alJ 1 1992k ICrowfher et alJ 
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120031) . This preliminary optical identifica- 
tion has been confirmed by a detailed study 
with the Chandra and Hubble observations 
(iBauer & Brandt! l2004h Furt her Chandra ob- 
servations dWang et al.1 12005) have confirmed 
the X-ray variability of IC10 X-l: the X-ray 
flux has been monitored over a period of half 
a day and it varied systematically by a factor 
of four. This suggested that IC10 X-l may be 
eclipsing accretion powered binary. An anal- 
ysis of a long time X-ray flux using Chandra 
data together if the X-ray monitor on board 
of Swift sa tellite revealed the o rbital period of 
32 hours dPrestwich et al.ll2~007l) . In this paper 
the authors analyzed the spectroscopic obser- 
vations of IC10 X-l at two different epochs 
and obtained a lower limit on the radial veloc- 
ity amplitude. This lead to a preliminary esti- 
mate of the mass of the BH in the system. The 
mass of the WR star was estimated to be at 
least 17 M . For the inclination of 90 degrees 
this implied a BH mass of 23 M . For smaller 
inclinations and higher masses of the WR star 
the estimate mass of the black hole was larger. 
Thus, IC10 X-l hosts one of the most massive 
stellar BHs known. A few months later a de- 
tailed spectroscopic study of this binary lead 
to me asurement of a detailed radial velocity 
curve dSilverman & Filippenkdl2b08l) . and the 
the analy sis in this paper confirm s the mass es- 
timate of lPrestwich et alj (120071) . 

Thus IC10 X-l is a binary consisting of a 
BH with the mass of at least 23 M accreting 
from a WR star with mass of 17M or more. 
If the mass of the WR is larger than the BH 
is also more massive. For the upper limit on 
the WR star mass of 35 M and the smallest 
possible inclination still allowing for eclipses 
the estimate of the mass of the BH is « 34 M . 

3.2. Future evolution 

What will be the future evolution of IC10 
X-l and what are its implications? The WR 
donor star in the binary is losing mass at 
through its strong wind. The mass loss rates 
from WR stars are highly uncertain. They 
have been estimated with various observation 
techniques dNieuwenhuiizen & de Jagerll 19901: 
Hamann & Koesterke 2000; Nugis & Lamers 



2000). We have estimated the expected mass 
loss rates from stars in low metallicity envi- 
ronment like the one in IC10 and calculated 
the final remnant mass as a fu nction of the ini- 
tial star mass, for details see Belczvnski et al. 
(2009). One clue to the mass loss rates comes 
from the existence of the BH with the mass of 
at least 23 M . In order to form such a massive 
BH the initial star must have been very mas- 
sive and the winds strength must not have been 
strong. In fact only the models where the winds 
a scaled down by at least 50% from the fiducial 
values REF lead to formation of BHs with the 
mass exceeding 23 M . 

Assuming that the current mass loss rate is 
not stronger than the one offered for the pro- 
genitor of the BH in the system we find that 
the WR star will lose a few solar masses be- 
fore it collapses and forms a BH. The mass of 
the BH formed in this collapse will not differ 
much from the mass of the collapsing star and 
will likely be * 13 M Q , assuming the present 
mass of 17 M . If the mass of the WR star is 
larger than the BH will be more massive. 

The typical lifetime of a 17 M WR star is 
several hundred thousand years, certainly not 
exceeding 10 6 yr. What is remarkable about this 
binary that it is a system that defies the "bot- 
tleneck" problem outlined in section 2. The 
current mass of the first formed BH is large 
enough to make the mass transfer from the 
companion stable. Thus the apparent problem 
in formation of BBHs outlined above is au- 
tomatically solved. The system is in a stable 
mass transfer state after the formation of the 
first black hole and its companion is likely to 
form another BH within the next few hundred 
thousand years. 

Also the binary should not be disrupted 
in the formation of the second BH. First the 
system is not going to loose a lot of mass 
when the BH is formed. Second, the kick ve- 
locity may be imparted on the newly formed 
BH. However, numerical estimates of such 
kicks show that they are smaller than in the 
case of neutron stars, and may reach at most 
150km s~ l . The relative orbital velocity of the 
components of the binary has been measured 
to be « 800km s -1 . The kick velocity would 
have to be larger than that in oder to disrupt 
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the system. Thus the system is not going to be 
disrupted by the formation of the second BH. 

Therefore in a few hundred thousand years 
IC 10 X-l will evolve into BBH system with the 
masses of at least 23 M and 13 M . The orbit 
will most likely be similar to the current one, 
however it may be altered by the natal kick. 
If the orbit remains the same then the merger 
time of the BBH binary will be w 3Gyrs. A de- 
tailed simulation with the kick velocity drawn 
from Maxwellian distribution with the width 
of 150km s _1 , has shown that only in 3% of 
cases the merger time of the BBH exceeds the 
Hubble time. 

4. The binary black hole population 

4.1. Estimate of the BBH merger rate 
density 

The estimate of the cosmic merger rate den- 
sity based on a single object may seem like fu- 
tile task yet in the field of gravitational wave 
astronomy this has already been done. The ob- 
servational estimate of the BNS merger rate de- 
pends crucially on the observations of a single 
object, namely the binary pulsar J0737-3039 
dBurgay et al]l2003h . 

In order to estimate the rate density we 
need to evaluate the volume in which IC10 X- 
1 could have been found, as well as the time 
in which it is observable. In order to estimate 
the volume we need to analyze critically the 
crucial observations that led us to the identifi- 
cation of the nature of IC 10 X-l, and find the 
most constraining one. IC10 X-l was identi- 
fied as a variable X-ray source, and X-ray ob- 
servations allowed to determine its orbital pe- 
riod. The second crucial observation was the 
spectroscopic measurement of the the radial 
velocity curve. This is the most constraining 
one: the spectroscopy of star is possible down 
to the apparent magnitude of my « 21. For a 
WR with the absolute magnitude of M » -5 
this corresponds to the distance modulus of 
25, i.e. the distance of R obs = 2Mpc. Thus 
the volume in which IC10 X-l is observable 
is V obs = 4nR 3 obs /3 ~ 33.5 Mpc 3 . In the fol- 
lowing we will assume conservatively that the 
entire sky has been surveyed for IC10 X-l like 



objects. Since IC10 X-l has been discovered 
and identified in X-rays, the duration of the X- 
ray bright phase defined the time it is observ- 
able. For IC10 X-l we assume conservatively 
that the observability time is equal to the life- 
time of the WR star, namely t obs = 0.5Myr. 
Therefore the local formation rate of IC 10 X-l 
like objects is 

<R = {Vt obs y l = 0.06Mpc- 3 Myr _1 . (1) 

Since nearly all of the systems that form will 
merge within the Hubble time this is also the 
estimate of the merger rate density. 

4.2. The expected BBH merger rate 

The gravitational wave signal from a coa- 
lescing binary in the inspiral phase depends 
on a single parameter, the chirp mass M - 
( m i m 2 ) 3 ^ 5 ( »z i + tn%f x ' , where m, are the 
masses of the components of the binary. For 
a given signal to noise ratio the range up to 
which a binary can be detected is D <x M 5 ^, 
and thus the volume where binaries can be de- 
tected scales as Vd et M 5 ^ 2 . The chirp mass 
for the BBH that will form from IC 10 X-l is 
« 14 M , while the chirp mass for the double 
neutron star (DNS) binary consisting of two 
1.4 M star is « 1.2M . Thus the the BBH 
that forms out of IC10 X-l from a distance 7.8 
times larger than the fiducial DNS binary. For 
the already finished S5 run of the LIGO de- 
tector the detectability range for DNS system 
was as 18 Mpc, so the detectability distance for 
a BBH with the chirp mass of m 14 M was 
Dbbh ~ 140 Mpc. We can now estimate the ex- 
pected BBH coalescences rate, as the volume 
observed was V B bh ~ 11.5 x 10 6 Mpc 3 . Thus 
the expected detection rate is 

N = 0.69Dl m M 5 ^ H R-lr o L (2) 

where D DNS = D^/ 18 Mpc, Mbbh = 
Mbbh 1 143 M q , R obs = R obs /2 Mpc, and t obs = 
t obs /0.5 Myr. One has to remember that this es- 
timate is based on a single object so in order to 
be fair we have presented explicitly the scal- 
ings. However we will argue below that this 
estimate of the detection rate should be consid- 
ered as a conservative lower limit. For a more 
detailed discussion see Buli k etal] (120081) . 
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5. Conclusions 

The detection rate obtained above has to be 
treated with some caution as it was obtained 
with the use of just a single object. However, 
we must stress that the assumptions made in 
calculating it are rather conservative. The life- 
time of the WR star has been assumed to 
0.5 Myr, yet for the higher masses of the star 
this time is even smaller, which makes the rate 
density go up. Also we have been quite conser- 
vative in determining the detectability distance 
for IC10 X-l like binaries. The spectroscopy is 
difficult even for a WR star at the distance of 
to IC10, but with the best telescopes available 
it might be possible up to 2 Mpc. If this dis- 
tance is actually smaller than the value we have 
assumed than the estimate of the rate density 
will also go up. We have also made an assump- 
tion that the entire sky has been surveyed for 
IC10 X-l like objects. This is also conserva- 
tive because if the actual volume surveyed was 
smaller than the rate density has to go up again. 
We have also neglected the upward correction 
on the rate due to the fact that IC10 X-l is an 
eclipsing binary. 

In the calculation of the detection rate we 
made a crucial assumption that the rate density 
estimates locally can be extrapolated to the dis- 
tance of hundreds of Mpc. This is an assump- 
tion that is also made in case of estimates of the 
DNS coalescence detection rate. On the scales 
of hundreds of Mpc the mean galaxy density 
may be a little lower than in the Local Group 
within the distance of 2Mpc, yet the effect can- 
not be very significant. It may decrease the de- 
tection rate by a factor of a few at most. The ex- 
pected detection rate depends strongly on the 
estimated chirp mass of the future BBH bi- 
nary. This in turn depends on the strength of 
the winds and the amount of wind mass loss 
from the WR before it forms a BH. However, 
we have an estimate of the wind from the ex- 
istence of the massive BH in the system, and 
also form observations. These two estimates 
imply that the mass of the compact object to be 
formed from the WR star should not be small. 
It should be stressed at this point that the un- 
certainty in the mass of the BH and the ensu- 
ing uncertainty in the chirp mass of the BBH 



is the most significant source of uncertainty in 
the estimate of the BBH coalescence rate. 

The astrophysical significance of this rate 
is quite interesting. First, it is much higher than 
the double neutron star coalescence rate. The 
rate is so large that the probability of detec- 
tion in the LIGO S5 and VIRGO VSR1 data 
is significant. If the data analysis of these ob- 
servations does not lead to a detection then the 
upper limits will certainly be very intriguing. 
The LIGO and VIRGO detectors are operat- 
ing now with the enhanced sesitivity. Thus, it 
is likely that the current observational runs of 
LIGO and VIRGO will finally be awarded with 
a detection and the first source of gravitational 
waves will be a coalescing BBH. 
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